Probability-preserving evolution in a non-Hermitian two-band model 
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A non-Hermitian PT-symmetric system can have full real spectrum but does not ensure probabil- 
ity preserving time evolution, in contrast to that of a Hermitian system. We present a non-Hermitian 
two-band model, which is comprised of dimerized hopping terms and staggered imaginary on-site 
potentials, and study the dynamics in the exact PT-symmetric phase based on the exact solution. 
It is shown that an initial state, which does not involve two equal-momentum-vector eigenstates in 
different bands, obeys perfectly probability-preserving time evolution in terms of the Dirac inner 
product. Beyond this constriction, the quasi-Hermitian dynamical behaviors, such as non-spreading 
propagation and fractional revival of a Gaussian wave packet, are also observed. 

PACS numbers: ll.30.Er, 03.65.-w, 03.75.-b 



I. INTRODUCTION 

Hermiticity of the Hamiltonian as the fundamental 
postulate in quantum mechanics guarantees the real 
eigenvalues and the conservation of probability. However, 
the recent discovery of Bender and Boettcher showed that 
Hermiticity of the Hamiltonian is not essential for a real 
spectrum [Ij. It has been proved that a non-Hermitian 
■pT-symmetric Hamiltonian can have real spectrum [2- 
|4|. Based on a time-independent inner product with a 
positive-definite norm, a new class of complex quantum 
theories having positive probabilities and unitary time 
evolution is established. The Hermitian and the non- 
Hermitian Hamiltonians seem to describe two parallel 
worlds, and much effort has been devoted to the con- 
nection between them |l|-[l6|. 

One of the ways of connecting a pseudo-Hermitian 
Hamiltonian with its Hermitian counterpart is the 
metric-operator theory outlined in y, providing a map- 
ping between them. However, the obtained equiva- 
lent Hermitian Hamiltonian is usually quite complicated 
[J, [l7|- Alternative ways, such as the interpretation of 
the non-Hermitian systems in the frameworks of scatter- 
ing and quantum phase transition, have been investigated 

[liil. 

In this work, we investigate the dynamics of a VT- 
symmetric pseudo-Hermitian Hamiltonian in the context 
of unbroken VT symmetry. We consider an exactly solv- 
able non-Hermitian VT model. It is a two-band tight- 
binding ring, with the non-Hermiticity arising from stag- 
gered imaginary potentials. It has been shown that such 
potentials can be realized in the realm of optics through 
a judicious inclusion of index guiding and gain/loss re- 
gions |23l - [26| . Recently, it was reported that the most 
salient character of the pseudo-Hermitian Hamiltonian, 
which is the VT symmetry breaking, was observed exper- 
imentally [23, [231 ■ Nevertheless, the reality of the spec- 
trum is not the unique common feature for the pseudo- 
Hermitian and the Hermitian systems in some cases. In 
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Ref. [231 it is pointed out that some non-Hermitian scat- 
tering centers, which consist of two Hermitian clusters 
with anti-Hermitian couplings between them, can act as 
Hermitian scattering centers, i.e. the S-matrix is unitary, 
or the Dirac probability current is conserved. The goal 
of the present work is to show the dynamical similarity 
between a non-Hermitian system and a Hermitian one. 
Intuitively, closely localized gain and loss potentials may 
be balanced with each other, or equivalently, the tempo- 
ral and spatial large-scale dynamics should be probability 
preserving. The Dirac inner product can be measured in 
an universal manner in experiments, hence it is of central 
importance to most practical physical problems. In this 
work we aim at investigating the dynamical behavior in 
terms of the Dirac inner product. Within the unbroken 
T^T-symmetric region, the eigenfunctions with different k 
are orthogonal spontaneously in terms of the Dirac inner 
product. This feature ensures the probability-preserving 
evolution of a state, which involves only one or two sub- 
bands with different k. In this sense, the non-Hermitian 
Hamiltonian acts as a Hermitian one without employ- 
ing the biorthogonal inner product. We also provide 
some illustrative simulations to show the occurrence of 
the fractional revivals and the slowly spreading of a wave 
packet. It shows that for certain special models, the non- 
Hermitian and Hermitian Hamiltonians can describe the 
same physics within a certain energy range. 

This paper is organized as follows. In Sec. [Hi we 
present a non-Hermitian T^T-symmetric two-band model 
and its exact solution. In Sec. lIIIi we investigate the Her- 
mitian counterpart of this model. In Sec. IIVI we show 
the quasi-canonical commutation relations and the quasi- 
Hermitian dynamics. In Sec. |Vl we demonstrate the re- 
sults for the system approaching to the exceptional point. 
Sec. lVIl is the summary and discussion. 



II. MODEL AND SOLUTIONS 

We consider a two-band model described by a non- 
Hermitian Hamiltonian H. It is a tight-binding ring 
with the Peierls distortions between nearest-neighboring 
sites and the additional staggered imaginary on-site po- 
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FIG. 1. (Color online) Schematic illustration of the Peierls 
tight-binging ring with staggered imaginary potentials (solid 
and empty circles). It is P7"-symmetric with respect to the 
OO' axis and invariant under the translation. 



tentials, which can be written as follows 
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where a\ is the creation operator of a boson (or a fermion) 
at the lih. site, with the periodic boundary condition 
ci2N+i = di- The hopping strengths, the distortion factor 
and the alternating imaginary potential magnitude are 
denoted by J, 5 and 7 (7 > 0), respectively. A sketch of 
the lattice is shown in Fig. [TJ In the absence of the stag- 
gered potentials or the Peierls distortion (with real poten- 
tials) , it is a standard two-band model and is employed to 
be a gapped data bus for quantum state transfer 30 32] . 
It is a T'T-symmetric model with respect to an arbitrary 
diameter axis. Here, without loss of generality, we define 
the action of time reversal and parity in such a ring sys- 
tem as follows. While the time reversal operation T is 
such that TiT — ~i, the effect of the parity is such that 
VajV ~ '^2N+i-i- ^PPlyiiiS operators V and T on the 
Hamiltonian ©, one has [T, H] =/= a.nd [P , H] ^ , but 
[VT, H] — 0. According to the non-Hermitian quantum 
theory, such a Hamiltonian may have fully real spectrum 
when appropriate parameters are taken. In the follow- 
ing, we will diagonalize this Hamiltonian and show that 
it has fully real spectrum. 

Beyond the VT symmetry, H is invariant under the 
translational transformation. Then taking the Fourier 
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where k — 2TTn/N, n £ [0, A^ — 1] is the momentum, the 
original Hamiltonian can be expressed as 



H = Y,H, 



(3) 



with 



(l-rf+(l + d)e^*'=)4Bfe+H.c 
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Hk 



Here Al and Bl are two kinds of creation operators of 
bosons (or fermions), resulting [Hk,Hk'] = 0. The op- 
erator Hk is non-Hermitian and can be readily written 
as 



Hk = -Cfe {akak - Pkh) 

by applying the linear transformation 

ttfc = ^fcAfe -f J^fcSfc, 
Pk = —VkAk + jlkBk, 

and 

ttfc = fikAl + VkBl, 
Pk ^-VkAl+ HkBl, 

where the spectrum is given by 
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where 6k and (pk are 
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with (j)k e [0, 27r], 5 > and Afe = 7/efe. 



The non-Hermitian operator Hk in Eq. ((4]) is in diago- 
nal form, since a/., a^, /3fc, and jSk are canonical conjugate 
operators, obeying the canonical commutation relations 



[ak: ak']± = [I3k: (3k']± = 0, 
[dk:ak']j- = [/3fc,/3fc']^ = 0, 



(9) 
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Therefore, the original Hamiltonian ([T]) is diagonalized. 
The method employed here is similar to that for the 
Hermitian two-band models [31"'33]. Nevertheless, the 
transformation in Eqs. ([5]) and ([6]) is no longer unitary 
under the Dirac inner product, since the canonical con- 
jugate pairs appearing in Eq. 1^ are not simply defined 
by the Hermitian conjugate operation, i.e. ak y^ aj. and 
/3/j 7^/3^, which is crucial in this work. 

We note that the spectrum ek consists of two branches 
separated by an energy gap 
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Obviously, it displays a full real spectrum within the re- 
gion of AJ^S^ > 7^. Beyond this region, the imaginary 
eigenvalues appears and the VT symmetry of the corre- 
sponding eigenfunction is broken simultaneously accord- 
ing to the non-Hermitian quantum theory. Interestingly, 
it occurs independently on the size of the lattice. No- 
tice that, when the onset of the VT symmetry breaking 
begins, the band gap vanishes, which is similar to that 
in a Hermitian two-band model. However, the dimeriza- 
tion still exists {S ^ 0), when the gap vanishes in this 
non-Hermitian model. In the next section, the further 
relationship between a non-Hermitian and a Hermitian 
two band models will be discussed. 



III. HERMITIAN COUNTERPART 

In this section, we would like to construct the equiva- 
lent Hermitian counterpart of the non-Hermitian model 
Eq. ([T|), which is a typical topic in the non-Hermitian 
quantum theory. In general, this can be done in the 
framework of metric-operator theory [3,|3]- Nevertheless, 
for the present model one can achieve this goal in a more 
direct way. This is due to the fact that the spectrum tk 
has an evident physical meaning. To demonstrate this 
point, we consider the model of a Peierls distorted tight- 
binding ring with staggered real potentials. The Hamil- 
tonian can be written as 

2^ 



+v,Y.^-\ih\h, (11) 



where h[ is the creation operator of a boson (or a fermion) 
at the /th site, with the periodic boundary condition 
&2JV+1 = &i- This Hamiltonian can be viewed as the Her- 
mitian counterpart, which will be shown in the following. 
By the similar procedure, taking the unitary transforma- 
tion 
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and He can be written in the diagonal form 



where the spectrum is 



efc = 2JeW(l~42)cos2 
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It is different from the situation of a non-Hermitian 
model, the coefficients C,k and ^k satisfy 



M' + M'^i. 



(15) 



which ensures the unitarity of the transformation in 
Eq. (J12p and the canonical commutation relation 
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Comparing two spectra e^ and e/c, one can see that they 
can be identical under the condition 



S^ ~ (7/2 J)' _ (52 + {Vj2J,f 
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Therefore, Hamiltonian He can be regarded as an equiv- 
alent Hermitian Hamiltonian of H. 

To illustrate this point, we consider a simple case of 
He with no energy gap A = and 7 = 7c = 2 JJ. Then 
the corresponding equivalent Hermitian Hamiltonian has 
the form 



he = - Jr 
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b%+i 



H.C. 



(18) 



which represents a uniform ring system with hopping am- 
plitude Je = JVl — S^. In Sec. |IV] we will investigate 
the wave-packet dynamics. It is noted that, although the 
spectrum for the non-Hermitian model is equivalent to 
that of a uniform ring, the distortions 6 and the imagi- 
nary potentials 7 are still nonzero and affect the dynam- 
ics in a balanced manner. 

We would like to point out that the method employed 
in this work is not universal as it depends on the obtained 
spectrum. We believe that the equivalent Hamiltonian 
He can be obtained by the standard metric-operator the- 
ory d, U. Actually, both methods have been used to 
another non-Hermitian model in a previous work [34| . 



IV. QUASI ORTHOGONALITY AND 
HERMITIAN DYNAMICS 



It is well known that the eigenstates of a non-Hermitian 
Hamiltonian can construct a set of biorthogonal bases 
in associate with the eigenstates of its Hermitian conju- 
gate. For the present Hamiltonian in Eq. ((IJ, eigenstates 
{ttk \0)Jk |0)} of H and eigenstates {al |0) ,/3^ |0)} of 
H^ are the biorthogonal bases of the single-particle in- 
variant subspace. This can be extended to the many- 
particle invariant subspace due to the canonical com- 
mutation relations in Eq. ^. On the other hand, the 
eigenstates of a non-Hermitian Hamiltonian are not or- 
thogonal under the Dirac inner product in the general 
case. However, we note that the eigenstates of the present 
Hamiltonian ([1]) are the eigenstates of momentum simul- 
taneously, which should lead to the orthogonality be- 
tween the eigenstates with different k in the Dirac in- 
ner product. This property is reflected by the following 
quasi-canonical commutation relations 
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Here the term "quasi" is the manifestation of the non- 
Hermitian nature of H in Eq. ©, which is represented 
in the absence of orthogonality between the eigenmodes 
of ctk and Pk- On the other hand, the rest "canonical 
commutation relations" makes the non-Hermitian system 
appear Hermitian to some extent. Similar relations and 
corresponding dynamical phenomena in a 'PT'-symmetric 
ladder system were presented in a previous work 



Now we turn to investigate the dynamics of such two- 
band model. Owing to the non-Hermiticity of the Hamil- 
tonian, the time evolution operator U (t) = exp {—iHt) 
is not unitary. To clarify the feature of the dynamics, we 
consider the time evolution of an arbitrary state. For the 
given initial state 

l'/'(0)>=^(/fcafc+3fc^fc)|0), (20) 

k 

we have 

= Y. (e^'^Vfettfe + e-'^'-'gkh) |0) . (21) 

k 

There are two types of probability, Pd (i) and Pb (i), 
in terms of the Dirac and biorthogonal inner product. 
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and 
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denote the Dirac and biorthogo- 



nal norms of the state 1-0), respectively. From the com- 
mutation relations Eq. ©, we have Pb [t] = 1, which 
is the aim of the introduction of the biorthogonal inner 
product. In contrast, Pp (i) is not unity and probably 
huge in some cases |34| . 

From the quasi-canonical commutation relations of 
Eq. p^ . we have 



1 + A^ 



^d(0 = E(IM' + I5'cI 

A; 

-f2^Afc|5fe/fc|sin(2efct-t-^fe), (24) 

k 

where A^ = 7/efc and ifk is a time-independent phase de- 
fined as e"^*-' = glfk/ \9kfk\- Obviously, the first term 
is time-independent while the second term represents 
a summation of periodic sinusoidal functions with fre- 
quency 2ek- In case of gkfk = (for each eigenmode 
k, the initial state does not comprise components of ak 
and /3fc simultaneously) and Xk being finite (the initial 
state does not comprise the component of e/j = 0, when 
the Hamiltonian becomes a Jordan block operator), the 
probability-preserving time evolution occurs. Neverthe- 
less, even in the case of gkfk 7^ 0, if Aj, ^ 1, the proba- 
bility slightly fluctuates around a certain constant, with 
the time evolution being quasi-probability-preserving. 



V. WAVE PACKET DYNAMICS 

Now we apply the obtained results to a more concrete 
case and then demonstrate the dynamic property of the 
system. We investigate the time evolution of the wave 
packet in the system with zero band gap. As mentioned 
above, it has been shown that the spectrum of the sys- 
tem is the same as that of a uniform ring, which can be 
regarded as the equivalent Hermitian counterpart. 

As an application of the obtained result, we consider 
the time evolution of a Gaussian wave packet (GWP) 
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with the central momentum k^ G [— ti", tt], centered at the 
A^Ath site, where |^) = a] |0) and J^i is the normaliza- 
tion factor. By using the inverse transformation from 






FIG. 2. The illustration of the time evolution of a Gaussian wave packet (solid line) with a = 0.1 and (a) fco = (b) fco = 3tt/8 
(c) fco — 7r/2 in a ring of TV = 100, S — 0.1 and 7 = 0.2 — 10~* ~ 7c (in units of J), where 7c — 2JS. We take t in units of Trcv 
from Eq. (|29|) in (a), and Tdv from Eq. (|31|l in (b) and (c). For comparison we also plot the same wave packet (hallow triangle), 
which evolves in a uniform ring of he from Eq. (|18p . One can see that the wave packet of fco = splits into several sub-GWPs, 
which almost have the same shape as the initial one and are referred as the fractional revivals. And those of fco — 3n/8 and 
7r/2 translate smoothly in the ring, where the latter behaves the non-spreading propagation. These figures show that the time 
evolution of the GWPs under the non-Hermitian H gives the quasi-Hermitian dynamical behaviors, which is similar to that 
under he. 
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FIG. 3. Plots of the square of the Dirac norm (solid line) from Eq. (|24|l in Fig. [2l For comparison we also plot the same wave 
packet (hallow triangle) in a uniform ring of he with the hopping amplitude of Je. We take t in units of Trev from Eq. (|29p in 
(a), and Tdi from Eq. (I3ip in (b) and (c). One can see that the Dirac norm fluctuates slightly and deviates little from unity. 
This is an obvious quasi-Hermitian dynamical behavior, which is in agreement with our predictions. 



the combination of Eqs. (I2|) and (O 
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the GWP of Eq. dlS]) has the form 
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with rj^ = — {v2TT-k) ■ It is a coherent superposition of 
eigenstates around k ~ 2fco in each band. However, in 
the case of |fco + tt/2\ :» 0, we have 
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?/fc/3fe|0), fco < -- or fco > 2 

Obviously, it satisfies the above mentioned probability- 
preserving condition of gkfk = 0, and then evolves as 
if in a uniform ring. On the contrary, in the case of 
fco ~ — 7r/2, we have |?7i!'/?7jr| ~ 1, i.e. two eigenmodes 
ak and f3k are both the main components of the state 
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FIG. 4. Plots of the square of the Dirac norm (sohd line) of 
the evolved Gaussian wave packet with a = 0.1, fco = — Stt/S 
and — 7r/2 in the same ring as Fig. [5] We take t in units of 
Tcir from Eq. pif) . For comparison we also plot the same 
wave packet (hallow triangle) as above in the uniform ring of 
he. One can see that the square of the Dirac norm fluctuates 
greatly and its average deviates from unity evidently, display- 
ing evident non-Hermitian behavior, which is in agreement 
with our predictions. 



FIG. 5. Plots of the square of the Dirac norm (solid line) of 
the evolved Gaussian wave packet with a = 0.1, fco = 7r/2 and 
— 7r/2 in the same ring as Fig. [21 but with 7 = 0.19 (in units of 
J). We take i in units of Tdr from Eq. 1)31^ . For comparison 
we also plot the same wave packet (hallow triangle) as above 
in a uniform ring of h^. One can see that the square of the 
Dirac norm fluctuates greatly and its average deviates from 
unity evidently, displaying obvious non-Hermitian behavior. 



simultaneously. From Eq. (|24p it is predicted that the 
dynamics of the wave packet should show extremely non- 
Hermitian behaviors. To demonstrate and confirm the 
analysis, we consider two typical cases of /cq = and 

7r/2. 

For fco = 0, at the instant i, we have 
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where Ti-ov is the characteristic revival time that can be 
estimated by [3a] as 



^ rnv — 



N' 



5^ 



2^2 

ttJc 



(29) 



It shows that the fractional revival occurs due to the 
approximate quadratic dispersion relation as if the wave 
packet evolves in the Hamiltonian h^. For fcg = ""/S, at 



the instant i, we have 
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which are the group velocity and circling period for a 
GWP of fco = 7r/2 in the effective ring. 

To demonstrate the above-mentioned results, numeri- 
cal simulations are performed. We plot the illustration 
and the square of the Dirac norm of the evolved GWPs of 
different /cq in this PT-symmetric two-band ring as well 
as a uniform ring for comparison in Fig. [51 We should 
notice that at the exceptional point, the gap disappears 
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FIG. 6. The illustration of the time evolution of the same GWPs on the same lattice as those in Fig. [2l but under open 
boundary conditions. We take t in units of Trcv from Eq. (|29|) in (a), and Tdr from Eq. (|31|) in (b) and (c). One can find the 
dynamics of the open boundary condition, which are the fractional revival and non-spreading propagation, to be similar as 
those in the periodic boundary conditions. 



and these two bands merge. Under this condition, the 
spectrum is the same as that of the effective uniform ring 
with the hopping amplitude being J^. For ko = 0, one 
can see that the GWP in the PT-symmetric ring has 
almost the same time evolution, which comes from the 
unequal distribution of the initial state on the two bands 
in the momentum space. For the wave packet with mo- 
mentum fco = 0; it mainly locates on the lower band of 
cxk around fc ^ and rarely locates on the upper band 
of /3fc, which satisfies the quasi-Hermitian condition of 
l5fc/fc| ~ 0. Under these circumstances, it can be treated 
as quasi-Hermitian and the dynamics of the wave packet 
is similar as well as in an effective uniform ring. And 
the Dirac probability of the wave packet slightly deviates 
from unity, which is plotted in Fig. [3] The situation is 
similar for a, ko — Stt/S wave packet, the Dirac proba- 
bility is also approximately conservative. For ko = 7r/2, 
although the wave packet consists of components from 
both bands of ak and /3fc, the quasi-Hermitian condition 
still fits. That is because, for the same eigenvector k, 
one component from the two different bands is almost 
zero and the other is finite while both zero on the bro- 
ken states of a,r and /3^. This meets the quasi-Hermitian 
condition and hence the specific GWP exhibits an analo- 
gous dynamical behavior as if in the effective Hermitian 
ring. 

We should notice that the Hermiticity of the evolu- 
tion on this T'T-symmetric ring depends on not only 
the Hamiltonian, but also the distribution of the wave 
packet on the two bands. At the exceptional point, only 
two eigenstates are broken. When the components of the 
GWP consist of neither the two states, the Dirac norm 
will probably be quasi-Hermitian. The numerical sim- 
ulations are plotted in Figures [5] and [31 It shows that 



the time evolution of fco = a-nd Stt/S for the unbroken 
Hamiltonian are about the same as those for the Hamil- 
tonian near the exceptional point. When the central mo- 
mentum kg changes, the distribution on the two bands 
changes (as plotted in Fig. 2]). The quasi-Hermitian con- 
dition of \gkfk\ ~ is invalid, then the Dirac norm de- 
viates from unity apparently. In an unbroken ring with 
7 — 0.19, for the wave packet of kg = tt/2 and —tt/2, 
which contains the two unbroken eigenstates of a^r and 
/3^ simultaneously, the quasi-Hermitian condition is no 
more satisfied and the wave packet behaves in a non- 
Hermitian way as plotted in Fig. [5l 



VI. SUMMARY AND DISCUSSION 

We have proposed a non-Hermitian T'T-symmetric 
two-band model, which consists of dimerized hopping 
terms and staggered imaginary on-site potentials. We 
have shown that such a model can have real spectrum 
and exhibit Hermitian dynamical behavior, obeying per- 
fectly probability-preserving time evolution in terms of 
the Dirac inner product. This fact indicates that the 
balanced gain and loss in a non-Hermitian system can 
result in quasi-Hermiticity. Apparently, such a dynam- 
ical behavior arises from the quasi-canonical commuta- 
tion relations in Eq. (jO]) . The essence is the translational 
symmetry of the model, which ensures the gain and loss 
to distribute homogeneously. It is presumable that simi- 
lar phenomenon occur in a two-band chain system. It is 
more difficult to get the analytical result when the open 
boundary is applied, compared to the periodic bound- 
ary condition. In this case, numerical simulations have 
be performed to compute the time evolution of a wave 



packet by direct diagonalization of the Hamiltonian. We 
plot the numerical result for the evolution of the same 
wave packet on the open chain in Fig. |6l It shows that 
the open boundary condition docs not affect the obtained 
result so much. Since an open chain is much more feasi- 
ble to realize in practice compared to the ring, our results 
can give a good prediction for the matter-wave dynamics 
in experiments. The recent observation of the breaking 
of VT symmetry in coupled optical waveguides [23, [2^ 
may pave the way to demonstrate the result presented in 



this work. 



ACKNOWLEDGMENTS 



We acknowledge the support of the National Basic Re- 
search Program (973 Program) of China under Grant No. 
2012CB921900. 



[1] C. M. Bender and S. Boettcher, Phys. Rev. Lett. 80, [19 
5243 (1998). [20 

[2] A. Mostafazadeh, J. Math. Phys. 43, 205 (2002). 
[3] A. Mostafazadeh, J. Phys. A 36, 7081 (2003). [21 

[4] A. Mostafazadeh, J. Phys. A 37, 11645 (2004). [22 

[5] Z. Ahmed, Phys. Lett. A 282, 343 (2001); 286, 30 (2001); 

64, 042716 (2001). [23 

[6] M. V. Berry, J. Phys. A 31, 3493 (1998); Czech. J. Phys. 

54, 1039 (2004). 
[7] W. D. Heiss, Phys. Rep. 242, 443 (1994); J. Phys. A: [24 

Math. Gen. 37, 2455 (2004). 
[8] H. F. Jones, J. Phys. A: Math. Gen. 38, 1741 (2005); [25 

Phys. Rev. D 76, 125003 (2007); Phys. Rev. D 78, 065032 
(2008). [26 

[9] J. G. Muga, J. P. Palao, B. Navarro, and I. L. Egusquiza, 

Phys. Rep. 395, 357 (2004). [27 

[10] C M. Bender, D. C Brody, and H. F. Jones, Phys. Rev. 

Lett. 89, 270401 (2002). 
[11] C M. Bender, Rep. Prog. Phys. 70, 947 (2007). [28 

[12] P. Dorey, C. Dunning, and R. Tateo, J. Phys. A: Math. 

Theor. 40, R205 (2007). 
[13] A. Mostafazadeh, Int. J. Geom. Meth. Mod. Phys. 7, [29 
1191 (2010). [30 

[14] M. Znojil, J. Phys. A: Math. Theor. 40, 13131 (2007); J. 

Phys. A: Math. Theor. 41, 292002 (2008); Phys. Rev. A [31 
82, 052113 (2010); J. Phys. A: Math. Theor. 44, 075302 
(2011). [32 

[15] C M. Bender, D. W. Hook, P. N. Meisinger, and Q. H. 

Wang, Phys. Rev. Lett. 104, 061601 (2010). [33 

[16] C M. Bender, D. W. Hook, P. N. Meisinger, and Q. H. 

Wang, Ann. Phys. 325, 2332-2362 (2010). [34 

[17] L. Jin and Z. Song, Phys. Rev. A 80, 052107 (2009). 
[18] L. Jin and Z. Song, Phys. Rev. A 81, 032109 (2010). [35 

[36 



L. Jin and Z. Song, Phys. Rev. A 83, 062118 (2011). 
L. Jin and Z. Song, J. Phys. A: Math. Theor. 44, 375304 
(2011). 

G. L. Giorgi, Phys. Rev. B 82, 052404 (2010). 
N. Moiseyev, Non-Hermitian Quantum Mechanics (Cam- 
bridge University Press, Cambridge England, 2011). 
K. G. Makris, R. El-Ganainy, D. N. Christodoulides, and 
Z. H. Musshmani, Phys. Rev. Lett. 100, 103904 (2008); 
Phys. Rev. A 81, 063807 (2010). 

Z. H. Musshmani, K. G. Makris, R. El-Ganainy, and D. 
N. Christodoulides, Phys. Rev. Lett. 100, 030402 (2008). 
R. El-Ganainy, K. G. Makris, D. N. Christodoulides, and 
Z. H. Musshmani, Opt. Lett. 32, 2632 (2007). 
A. Ruschhaupt, F. Delgado and J. G. Muga, J. Phys. A: 
Math. Gen. 38, L171-L176 (2005) 

A. Guo, G. J. Salamo, D. Duchesne, R. Morandotti, M. 
Volatier-Ravat, V. Aimez, G. A. Siviloglou, and D. N. 
Christodoulides, Phys. Rev. Lett. 103, 093902 (2009). 
C. E. Riiter, K. G. Makris, R. El-Ganainy, D. N. 
Christodoulides, M. Segev, and D. Kip, Nature Physics 
6, 192 (2010). 

L. Jin and Z. Song, Phys. Rev. A. 85, 012111 (2012). 
S. Yang, D. Z. Xu, Z. Song, and C. P. Sun, J. Chem. 
Phys. 132, 234501 (2010). 

M. X. Huo, Y. Li, Z. Song, and C. P. Sun, Europhys. 
Lett. 84, 30004 (2008). 

B. Chen and Z. Song, SCIENCE CHINA Physics, Me- 
chanics & Astronomy 53, pp. 1266-1270 (2010). 

W. P. Su, J. R. Schrieffer, and A. J. Heeger, Phys. Rev. 

B 22, 2099-2111 (1980). 

X. Z. Zhang, L. Jin, and Z. Song, Phys. Rev. A 85, 

012106 (2012). 

L. Jin and Z. Song, Phys. Rev. A 84, 042116 (2011). 

R. W. Robinett, Phys. Rep. 392, pp. 1-119 (2004). 



